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Abstract 

In layered polar semiconductor BiTel, giant Rashba-type spin-split band dispersions show up 
due to the crystal structure asymmetry and the strong spin-orbit interaction. Here we investigate 
the 3-dimensional (3D) bulk band structures of BiTel using the bulk-sensitive /^-dependent soft 
x-ray angle resolved photoemission spectroscopy (SX-ARPES). The obtained band structure is 
shown to be well reproducible by the first-principles calculations, with huge spin splittings of ~300 
meV at the conduction-band-minimum and valence-band- maximum located in the k z = tt/c plane. 
It provides the first direct experimental evidence of the 3D Rashba-type spin splitting in a bulk 
compound. 
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Symmetry of the crystal lattice plays an essential role to determine the electronic struc- 
tures of solids. Among them, the spatial inversion / is an important symmetry since it 
connects the momenta k and —k, which are related also by the time-reversal symmetry T. 
In the latter case, the spin is also reversed, and hence the interplay of I- and T symmetries 
in the presence of the spin-orbit interaction (SOI) is an important issue. Investigation of SOI 
in noncentrosymmetric materials originates in Dresselhaus and Rashba effects jl, 2|, which 
evaluate the spin-split band structures in zincblende and wurtzite semiconductors. More 
recently, the Bychkov-Rashba model 3[ was proposed as a suitable explanation of the prop- 
erties of the 2-dimensional electron gas with perpendicular electric field E z applied. In this 
model, the electron with momentum k and spin a feels an effective magnetic field perpen- 
dicular to k and E z , as represented by the Hamiltonian H^r — A a ■ (E z x k), with coupling 
constant A. Such SOI effects in inversion asymmetric systems introduce various spin-related 
electric properties, such as spin Hall effect |4j and spin-galvanic current and con- 

sequently, occupy an important position in the spintronics field. Unconventional quantum 
electronic phase is also expected to appear in noncentrosymmetric superconductors, induced 
by the mixing of spin-singlet and spin-triplet states js, 9|. From the viewpoint of elucidating 
the SOI-induced novel phenomena and spintronic functions, the development of compounds 
with high-energy-scale spin-splitting is highly desirable. 

Pursuit for giant spin-splitting has been mainly focused on the Shockley state of single 
crystals composed of heavy elements until now 10, lU|, which is characterized by the clear 
structural asymmetry at the topmost surface cutoff. Indeed, the huge Rashba-type splitting 
was observed in Bi/Ag(lll) surface alloyed system 12j, where the fc-linear spin splitting 
similar to the Bychkov-Rashba dispersion E ± (k) = h 2 k 2 /2m* ±aj{\k\ (m*: effective mass of 
electron) had been demonstrated by angular- resolved photoemission spectroscopy (ARPES). 
Here, «r is the Rashba parameter indicating the strength of Rashba-type SOI, and «r = 3 
eVA is achieved in Bi/Ag(lll) system. Most recently, on the other hand, the comparable 
size of spin splitting («r = 3.8 eVA) has been found to be realized in a bulk polar semi- 



conductor, BiTel 



13]. Due to the peculiar trigonal polar structure (P3ml) arising from the 



alternative stacking of Bi- Te- and I- triangular networks [see Fig. 1(a)], it also shows a 
Rashba-type spin splitting similar to the ones observed in 2-dimensional surface Shockley 
states. Quantitative agreements between the (spin-resolved) ARPES and the first-principles 
calculations, together with the optical study on the bulk state, strongly suggest that the spin 
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splitting is derived from bulk BiTel {13I . [bll I15I ] . Nevertheless, the direct observation of 3- 
dimensional (3D) spin-splitting beyond Bychkov-Rashba model has not been investigated to 
date. Recent theoretical discovery of pressure-induced band inversion and topological tran- 



sition in BiTel 



161 ] further increases the importance of directly confirming the bulk band 



structure in detail. In this paper, we investigate the band structure of BiTel throughout the 
whole 3D Brillouin zone, by utilizing the photon-energy (hv) dependent bulk-sensitive soft 
x-ray ARPES (SX-ARPES). The result gives the first direct experimental evidence of 3D 
Rashba-type spin-splitting, which leads to a deeper understanding of SOI in solids. 



Single crystals of BiTel were grown by Bridgman method [13j. The electron-carrier den- 
sity was estimated to be 4.5 x 10 19 cm -3 from the Hall coefficient, showing the n-type de- 
generate semiconducting behavior, /izz-dependent SX-ARPES experiments were performed 
by using a Gammadata-Scienta SES2002 hemispherical electron analyzer at the undulator 
beam line BL17SU of SPring-8. The system adopts a grazing-incidence geometry (<5°) 
spectrometer that enables ARPES measurements to be highly efficient and allows an ac- 
curate determination of the k x — k y momentum due to negligible modification of the k z 



momentum 
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181 ] . The sample was cleaved in situ under an ultrahigh vacuum of 1x10 
Pa. All experiments were carried out at 50 K. The Fermi level of the sample was referred 
to that of the polycrystalline Au electrically connected to the sample. In order to avoid 
the contamination of Auger lines, soft x-ray between 530 eV and 590 eV were used. Both 
right and left circularly polarizations were used to avoid any polarization selection rule. The 
momentum region of /izz-dependent SX-ARPES measurement is shown in Fig. 1(b) together 
with the Brillouin zone of BiTel. All measurements were performed in the red translucent 
plane, which nearly corresponds to the T-M-L-A [k x — k z ) plane. Figure 1(c) shows the 
momentum cuts of measurements we can access by using the photon energies of 530-590 eV 
at intervals of 5 eV. The red solid (blue broken) curves correspond to the measurements in 
which the energy resolution was set to 100 meV (150 meV). 

SX-ARPES images along k x direction are shown in Figs. 2(a)-(c) for hv = 540, 565 
and 585 eV, which correspond to the k z values of 0, 0.567r/c and tt/c, respectively [see Fig. 
1(c)]. Peak positions of SX-ARPES intensity with error-bars estimated from momentum 
distribution curves (MDCs) or energy distribution curves (EDCs) are plotted as black-circle 
markers on SX-ARPES images, to show the dispersions of valence bands (VBs) below the 
band gap. At k x = 0, six VBs are observed below E-p, which mainly consist of Te and I 
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hp orbitals. Some of them further show splitting in k x 7^ region due to the SOI, as will 
be discussed later. Observed band structures show clear /^-dependence reflecting the 3D 
bulk electronic structure. Focusing on the highest valence band (HVB), its top lies at the 
binding energy (Eb) of 1.0 eV at k z = while it moves upward and forms the VB maximum 
(VBM) of Eb = 0.48 eV at k z = tt/c. The conduction band (CB) above the band gap, on 
the other hand, seems to be located above E-p at k z = 0, but gradually moves below the 
Fermi level and appears at around k x = for k z = 0.56tt/c and tt/c. The energy position of 
the conduction band minimum (CBM) is Eb ~ 0.22 eV at k z = tt/c. Thus in the bulk band 
structure of BiTel, both the VBM and CBM are located within k z = tt/c plane, forming the 
minimum gap of ~ 0.26 eV. 

To compare with the bulk band dispersions recorded by SX-ARPES, the results of the 



relativistic first-principles calculations 



19] for k z = 0, 0.567r/c and tt/c are shown in Figs. 



2(e)-(g). The energy axis of the calculated data is rigidly shifted to match the energy of VBM 
at Eb = 0.48 eV obtained by SX-ARPES. Calculated VBs consist of 12 bands because of 
spin splitting by SOI. All bands are spin degenerate along T-A line {k x = k y = 0) due to the 
time-reversal symmetry and the P3ml symmetry of the crystal. Although some deeper VBs 
at Eb = 2 ~ 6 eV tend to have somewhat lower binding energies compared to the SX-ARPES 
data, the shapes of the calculated bands coincide well with those obtained experimentally. 
The spin splittings for these deep VBs are also partly observed by SX-ARPES, e.g. the 
deepest VB showing the energy split of ~ 100 meV around k x = ±0.1 A~ x in Fig. 2(c). The 
calculated CB bottom (HVB top), on the other hand, gets deeper (shallower) in binding 
energy with increasing k z , and the minimum gap is formed at k z = tt/c. Simultaneously, 
the spin splittings of the CB and HVB become bigger, showing the huge band splitting of 
about 300 meV realized at CBM and VBM in the k z = tt/c plane. These are found to 
show the fc||-linear spin splitting in a small fcy region (k\\ = y/k% + k^), as similarly reported 



in wurtzite cases 



20|. Also in SX-ARPES image at k z = tt/c [Fig. 2(c)], the similar 



spin splitting of the HVB is clearly seen in the green rectangle, as will be discussed later 
quantitatively together with that of CB. 

In Fig. 2(d), the band dispersions along T-A line obtained by SX-ARPES using hv of 530 
~ 590 eV are shown, to elucidate the k z - dependence of the band structure. The calculation 
as shown in Fig. 2(h) accords well with the band structure obtained by SX-ARPES. The 
CB, HVB, and the deepest VB show moderately strong /^-dependence (0.4 ~ 0.8 eV) as 
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compared to other bands (< 0.2 eV), since they are mainly characterized by Bi 6p z , Te 5p z , 



and I 5p z , respectively 



141 ] . Still they are apparently weaker than the in-plane dispersions. 



This anisotropic electronic structure reflects the layered crystal structure of BiTeh which 
is also confirmed by the anisotropic effective mass of m z /m x ~5 by optical study 15j. On 
the other hand, since the CB crosses E-p between k z — and tt/c, it should show up as the 
n-type closed spin-split Fermi surface pockets centered at k z = ir/c, whose topology should 



be characterized by 3D-Rashba model 



2l|. 



To quantitatively compare the results of spin splitting obtained by SX-ARPES with that 
of the calculation, a close-up around the minimum gap, corresponding to the green-rectangle 
region in Fig. 2(c), is shown in Fig. 3(a). The peak positions estimated by EDCs [shown 
in Fig. 3(b)] and MDCs are represented as circles and squares for HVB, while the guides 
for eyes are shown as green broken curves to trace CB. ARPES intensity of the CB is much 
weaker than that of the HVB in Figs. 3(a) and 3(b), because of the small photoabsorption 
cross-section of the Bi 6p orbitals. The relativistic first-principles calculations as shown in 
Fig. 3(c) indicates the Rashba-type spin splitting of CB and HVB due to the SOI under the 
P3ml symmetry, which are very similar to the SX-ARPES results. To evaluate the strength 
of the Rashba-type spin splitting in analogy with the 2DEG Bychkov-Rashba model, we 
quantify the aforementioned Rashba parameter «r. Here we estimate «r by using the 
momentum offset at the CBM (/c^ BM ), the binding energy of the CBM (-Ecbm) and the outer 
Fermi momentum (k F nt ). Most simply, by assuming the parabolic band dispersions, they are 
related by a R = h 2 k% BM /m x and E C bm = h 2 (k F ut -k% BM ) 2 /2m x . Our SX-ARPES spectra at 
k z = n/c shows the "W" like shaped CB with \k$ BM \ = 0.04 ±0.01 A -1 , \k F ut \ = 0.10 ±0.01 
A" 1 and E CB m = 220±30 meV. Thus we can evaluate the Rashba parameter as «r = 4.9±1.6 
eVA. The value of lies within error bars of the value previously obtained from the CB 
in the surface accumulation layer of BiTel, and gives one of the strongest Rashba-type spin 
splitting as reported. Spin splittings can be also observed at the HVBs in Fig. 3(a). The 
offset momentum at the VBM is = 0.05 ±0.01 A -1 , whereas the spin splitting energy 

is estimated to be ~ 220 meV. These are fairly similar in size with the splitting of CB as 
mentioned above. 

This huge spin-splitting has been successfully explained by means of k ■ p perturbation 
theory [14J]. According to the theory, the SOI-induced band splitting near the Fermi level 
( i.e. CBM and VBM) is strongly affected by the unusual electronic structures around its 
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band gap; most importantly, the band-gap size and the same symmetry characters of CB 
and HVB due to the negative crystal field splitting of HVB. In BiTel, the VBM is mainly 
composed of Te 5p orbital whereas the CBM is characterized by Bi 6p. However, both the 
sets of bands share the same symmetry character within C% v space group, thereby leading to 
non- vanishing matrix elements, coupling them through the second-order perturbative k ■ p 
hamiltonian. The fairly narrow gap realized in BiTel (~0.36 eV from optical study [jjjl) 
combined with the large atomic SOI of Bi further enhances this coupling. These microscopic 
electronic structures should be also taken into account for designing the SOI effect in bulk 
compounds, besides the inversion-asymmetric crystal structure and the atomic spin-orbit 
coupling. 

Regarding the minimum band gap, it is estimated to be ~ 260 ± 50 meV at k z = n/c 
from the SX-ARPES data [Fig. 3(a)]. Coincidentally it is quite close to the calculated 
value (~ 0.28 eV), but considerably smaller compared to the result of optical spectroscopy 
showing the minimum band gap of ~ 0.36 eV 15[ . This discrepancy between the bulk may be 
arising from the modification of the electronic structure in the sub-surface (<2 nm) region. 
It can be qualitatively confirmed from the depth of the CBM. From the band calculation, 
E-p is expected to be around 142 meV above the CBM, if we assume that the Hall number 
hh = 4.5 xlO 19 cm" 3 correctly represents the bulk electron density. The present SX-ARPES 
data, nevertheless, indicates Eqbm ~ 220 meV, which is deeper than the calculation by 80 
meV. This difference may suggest that SX-ARPES data probes the 3D band dispersions 
from the boundary region between bulk and sub-surface (the probing depth of hv = 550 eV 
ARPES is typically ~ 1 nm), which is slightly shifted downward by the band bending effect 



131 ] . Recently, the formation of the charge accumulation layer and the resulting band-gap 



221 ] . Similar situation may be also realized in the 



narrowing have been reported in CdO 
sub-surface region of BiTel. 

Finally we would like to discuss the topology of the 3D Rashba-split bands. The calculated 
Fermi surface (Eb = 0) and the constant-energy-surface (CES) of HVB at Eb = 0.58 eV are 
shown in Fig. 4(b,c). E B = 0.58 eV corresponds to 100 meV below the VBM (E YB m ~ 0.48 
eV) and 70 meV above the band crossing point (Eb ~ 0.65 eV), as shown in Fig. 4(a). 
Since the ARPES intensity of the CB forming the Fermi surfaces is very weak, here we focus 
on the CES of HVB at E B = 0.58 [Fig. 4(c)]. This CES indicates a doughnut-like topology 
with the strong trigonal deformation, which is represented by a set of two separate distorted 
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ellipsoids of origin-symmetry, along T-A-L-M plane. This topology of CES reflects the 3D 
Rashba-type splitting of VB as a consequence of SOI in P3ml symmetry. To confirm this 
experimentally, we plot in Fig. 4(d) the /^-dependent SX-ARPES intensity at E B = 0.58 
eV with the energy window of ±4 meV. The upper region of the image (tt/c < k z < 2ir/c) is 
obtained by symmetrizing the lower region (0 < k z < tt/c), which corresponds to the data 
recorded using hu = 540 ~ 585 eV. The distribution of intensity, centered at A-point, is 
in a good correspondence with the cross-section view of calculated CES (indicated by red 
curves). It thus shows the direct evidence of the band topology with the 3D- Rashba-type 
splitting, as shown in Fig. 4(c), indeed realized in the bulk BiTel. Since the Fermi energy 
lies slightly above the band crossing point of the conduction band, the actual Fermi surfaces 
in the bulk should be consisting of two components, inner and outer ones as shown in Fig. 
4(b). In its cross-section view along T-A-L-M plane, two distorted ellipsoids slightly intersect 
with each other in the vicinity of A point. Thus, the magneto-transport properties of this 
compound should be dominated by the electrons forming such 3D spin-split inner and outer 
Fermi surfaces. The recent optical study on series of dopant-introduced BiTel shows that 
the chemical potential of this compound can be continuously varied from CBM to above 
the crossing point [ijj]. Together with this feasibility of carrier control, BiTel possesses the 
potential to realize various spin-dependent transport/optical properties {4] and spintronics 
as well as the new type of SOI-induced unconventional superconductivity 

In summary, we investigated the Rashba-type spin splitting in bulk of BiTel using /in- 
dependent SX-ARPES. Band dispersions indicated a moderate k z - dispersion with the min- 
imum band gap at k z = tt/c, showing a good agreement with the relativistic first-principles 
calculations. The conduction band minimum and the valence band maximum located within 
the k z = tt/c plane both showed a huge splitting, thus offering the direct evidence of giant 
Rashba-type spin splitting in bulk. The constant-energy surface at the valence band top 
indeed showed the doughnut-like topology, confirming the 3D Rashba-split Fermi surface 
realized in the bulk BiTel. These results show the potential of BiTel to realize the strong 
spintronic functions. 
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FIG. 1: (color online) (a) Crystal structure and (b) Brillouin zone of BiTel. (c) Momentum cuts of 
measurement in k x — k z (T-A-L-M) plane, corresponding to the photon energies of hv = 530 ~ 595 
eV (5 eV steps). Those recorded by energy-resolution of AE = 100,(150) meV are shown by solid 
(broken) curves. 

FIG. 2: (color online) (a)-(c) SX-ARPES images along k x direction by using hv = 540, 565, 
and 585 eV, respectively. Black circle markers with error bars indicate the peak positions of the 
SX-ARPES intensity. The green rectangle in (c) shows the CBM and VBM located around A 
point, (d) Out-of-plane band dispersions along T-A line obtained by SX-ARPES using hv = 530 
~ 590 eV, which correspond to k z = —0.22ir/c ~ l.ll7r/c. (e)-(g) Band dispersions along k x 
direction obtained by calculation corresponding to k z = 0, 0.567r/c, and tt/c, respectively, (h) 
Band dispersions along T-A line obtained by calculation. 

FIG. 3: (color online) (a) SX-ARPES image around the CBM and the VBM at k z = n/c (hv = 585 
eV). Open circle (square) markers on the HVB are the peak positions of the EDCs (MDCs), whereas 
the broken curves overlaid on the CB are guides for eyes. The corresponding EDCs are shown in 
(b). (c) The band calculation near the Fermi level at k z = tt/c along k x direction, showing the 
CBM and VBM. 

FIG. 4: (color online) Calculated band dispersion at k z = tt/c (a), 3D views of Fermi surface 
(Eg = 0) (b), and constant energy surface (CES) at Eb = 0.58 eV (c). (d) SX-ARPES intensity 
image at Eb = 0.58 eV in k x — k z plane, which are recorded by using hv = 540 ~ 585 eV. The 
right inset (the red curves in the main panel) indicates the calculated (cross-section) view of CES 
at E B = 0.58 eV. 
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